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Introduction 


All structural components of the solid rocket motor case of the space 
shuttle are considered to be fracture critical. It is also the present 
plan to reuse the solid rocket motor case for a designated number of 
missions. The expected number of missions and operations such as the 
tests on the case between the missions are accounted in the projected 
design life of the structure. A fracture control plan is necessary 
because fracture critical components are being reused. 

In particular, this report is concerned with the fracture control of 
the membrane of the six cylindrical segments that are considered to be the 
most critical of all structural components of the case. The developed 
procedure can, however, be used for all similar structures. During each 
mission, significant loads are applied to these six cylindrical segments 
during the flight and "slap down 11 operations. The applied stresses from 
all other events during the mission are considered not significant enough 
to result in cyclic or time dependent crack growth.* - If the test or analysis 
indicate the possibility of other critical loading events they can be in- 
cluded in the fracture control plan by extending the reported analysis. 
Before each mission, the cylindrical segments are also subjected to a 
proof test. The loads applied during the proof tests can result in sig- 
nificant amount of crack growth. As a preventive measure to reduce the 
effective depth of cracks, the thickness of the membrane is reduced by a 
selected amount between two missions. While the effective depth of crack 
is reduced, the operation has the effect of increasing the applied 


stresses. This necessitates a larger initial thickness of the membranes 
than that would be designed without this particular plan for fracture 
control. 

Therefore, any design of the membrane of the six cylindrical segments 
of the sd lid rocket motor case must arrive at an initial wallthickness "t" , 
the thickness 'A t f that will be decreased between each mission and the 
proof load factor 'K^ 1 . For example, a large value of initial wall thick- 
ness results in increased reliability, but results in the need for increased 
propellant, increased cost of operation and reduced pay load capability. 

On the other hand, a small initial wall thickness increases the probabi- 
lity of failure and the resulting loss of the shuttle vehicle and the pay 
ioad. Therefore, there is a need for optimizing the initial wall thick- 
ness. Similar arguments can be presented to explain the need for selecting 

the other design variables such as 'AL f and 'K 1 by optimizing the desired 

P 

objective function o£ cost and weight. 

In general, these design variables depend on the probability distribu- 
tion for the initial flaw sizes present in the membrane, applied stresses 
during the use of the vehicle, crack growth characteristics of the 
material, fracture control plans, specified reliability bounds^ weight 
and cost considerations. The report describes a reliability-based pro- 
cedure that can be used to select the design variables of SRM by using 
probabilistic fracture mechanics and cost or weight considerations. 

Method of Approach 

As discussed in reference 1, careful NDl' techniques can detect 

initial cracks greater than the surface length of c^ - 0.1 inch and surface 

2 



depth of a = 0.5 C with 1007* success. It has been claimed that cracks 

corresponding to surface length c^ = 0.1 inch can be identified 100% of 

the tine. If the corresponding maximum depth is 0.05 inch' thepe is no 

possibility of any initial cracks of depth larger than 0.05 inch. Such 

2 

an initial crack depth distribution can be analytically represented by 
Johnson distribution. The density function for the probabilistic 
model is written as follows 

+ (x T l)) _ • 

e < < e + A , rj > 0 , -<» < y < co , A>0 , (1) 

The four parameters of the distribution are j|, e anc * y 

This probability distribution for initial crack depth changes after 
each mission, proof test and the material removal from the wall thickness. 
The change in distribution after each mission and proof test is due to 
the crack growth resulting from the applied stresses. This crack growth 
also depends on the lengths of the crack that are already present and 
the material properties that are responsible for the crack growth. In 
this analysis, the applied stresses and material properties are assumed to 
be know- deterministically. If the initial crack length were also known 

deterministically the crack length after each use can be determined from 

3 4 5 

equations such as Paris’ equation , Foreman's equation or Collipriests 

equations. Because initial crack lengths are not known deterministically, 
crack length after each use of the vehicle is again another probabilis- 
tic distribution that has to be determined. 
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E <0 = 

a o 
o 




(a - e) (X - a + e) 


exp 
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The cumulative density function for crack length after 'n 1 uses is 
denoted by F(a n ). This represents the probability that a^ ^ A after n 
uses. Each use is defined as one flight, slap down, proof test and 
material removal. In this analysis "slap down" effects have not been 
considered. The crack growth due to "slap down" effects can be considered 
in a similar way. Also, crack growth due to time related effects such as 
stress corrosion have also been neglected. 

If F( a n ) is known the probability distribution for the stress in- 
tensity factor (K) can be obtained from the knowledge of the applied 
stresses. The probability distribution F(K q ) for stress intensity factor 
can be used to estimate the probability failure (P^) which is the pro- 
bability of stress intensity factor greater than or equal to the criti- 
cal stress intensoty factor during the projected design life of the 
structure. The critical stress intensity factor is denoted by K . In 
this analysis, stresses and the material properties are assumed to be 
known deterministically. However*, the applied stress changes after each 

use due to material removal. Therefore, the probability of failure can 

c c 

be expressed as the probability of a > a . In this expression a is 

n — 

the critical crack depth that can be obtained from the critical stress 
intensity factor and the applied stress. This relationship between the 
stress intensity and the applied stress is discussed in the next section* 

Stress Intensity Factor 

For the analysis of the stress intensity factor in the membrane, 
an infinite plate model with elliptical surface flaws that are oriented 
perpendicular to the applied stress has been assumed. The relationship 
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between the stress intensity, the applied tensile stress and crack depth 
is given by 


K 


where 



( 2 ) 


2 

Q (t) = * 2 - °- 212 2 ? (3) 

In this equation, crY is the yield stress and 0 is a function of 

?- 

the ratio of crack depth to crack length (a/c). Variation ft with (a/c) 
is given in reference 1. 

Because the crack depth (a) is a random variable the stress intensity 
factor K is also a random variable. In general, both crack depth a and 
crack length c are random variables and there is a need for a joint dis- 
tribution for a and c. In this analysis, only the crack depth is con- 
sidered as the random variable. It is also assumed that the probability 
distribution for crack depth ’a 1 is known initially and is given by a 

Johnson S, distribution. The density function for the distribution is 
b 

given in equation (1). This probability distribution for crack depth 
changes with use. The next step will be to determine the change and 
the new probability distribution after each flight and proof test. 


Probability Distributions for Crack Depth After Use 
The following symbols are used to properly account for the changes 
in probability distributions. 
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f(a^): Probability density function for the initial crack depth 

F(a^): Cumulative distribution function for initial crack depth 

F(a ): Cumulative distribution function for initial crack depth 

op 

after the first proof test 

F(a )• Cumulative distribution function after N flights and (N-fl) 
tests 

F(a ): Cumulative distribution function after N flights and N proof 

np 

proof tests. 

F(^n) : Cumulative distribution function after material removal from 

the wall thickness. 

Similarly, density functions are denoted by lower case f f f . As 
discussed before, 'slap down' effects are not considered in the analysis 
but can be included by following a similar procedure. 

The rate at which crack depth increases is assumed to be given by 
Paris' equation. Then 


(AK) n 

where C and n are empirical constants. Alternately, the rate of crack 
growth can be assumed to be given by Foreman's equation of Collipries t ' s 
equation if they are found to represent the situation more accurately. 
For example, Collipriest ' s equation can be written as follows: 


da 

exp 


An K - An AK 

c O . , 7 

n arc tan h \ 


+ In 


£n K 


c exp 


+ £n K \ 

i 


An AK - 


An K (1 - R) + An AK 
c o 


An K (1 - k) - An K 
c Q 


( 4 ) 


where n is an empirical constant. These equations can be used to obtain 
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crack depth after N+l uses if the crack depth after N uses and N proof 
tests are known deterministically, i.e., 

^+1 = Vt-l Wp* 

Similarly, crack depth after the proof test can be determined from equation 
(3) or (4) if the crack depth before the proof test is known determinis- 
tically, i.e.. 


a 


UP 


a fa } 
NP 1 N J 


( 5 ) 


These functions represented by equations (5) or (6) can be determined 

analytically or in the form of quadratures from equation (3) or (4). 

From equation (5), can* be obtained for every known value of a^. 

Similarly, a^ can be obtained for every known value of a^ from equation 

(6). However, both a and a are random variables in the present 

NP N 

analysis. In this case equation (5) can be used to obtain the proba- 
bility distribution for a n+ ^ if the probability distribution for a^ is 
known by using the principle of transformation of random variables. It 
should be noted that all equations similar to (5) or (6) involving crack 
depths are increasing functions. This property is useful in transforming 
the random variables. 

For example, the probability density function for a ^ can be 
written as follows 


f(a n+1 ) = f[a n+l{ a np >] | 


(7) 
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similarly 


da 


f (a ) = f [a {a } ] 
np np n 


da 


np 


( 8 ) 


Equations (7) and (8) can be written for every value of n from 
zero to the projected number of uses. 

Details of obtaining these equations for the membrane of the SRM 
with the expression for stress intensity given by equation (2) and Paris 1 
equation for crack growth is discussed in Appendix I. 

The next step is to obtain a tool for change of probability distri- 
bution due to the material removal from the wall thickness. 


Material Removal and the Change of Probability Distribution 
Due to material removal after each use the effective crack depth 
is reduced by ^ t 1 . Thus new crack depth is 

a = a -At 
n n 



It is assumed that is a constant. Thus, by- using the principles 

( 2 ) 

of transformation of random variables , the probability density func- 
tion for a. can be written as follows, 
n 


p(a ) = f ( a +4*0 
n n 


( 10 ) 


In this equation, p(a ) represents the density function for a and f 

n n 

represents the functional form of the probability density function for 


a . 
n 
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Probability of Failure 


By following the method discussed in the preceding two sections;, 
probability density function for crack depth can be obtained after every 
flight, proof test and material removal. From the density function, 
cumulative probabilities can be obtained by integration. Integration 
after the transformation of variables as discussed in equations (7), 

(8) and (10) needs the determination of appropriate limits of integra- 
tion consistent with the transformation of variables. This is also 
discussed in the Appendix I. If F(a^) represents the CDF after n flights 
and n proof tests the probability of failure is given by the probability 

of a > a . 
n c 

It is to be noted that the probability of failure changes with 
different selections of the initial wall thickness t, increased loading 
due to proof test, the material removed At and the number of designated 
missions. The increased loading due to proof tests is denoted by a factor 
Kp. A cost function or a weight function can be formulated from this 
knowledge of probability of failure and other related unit-cost or 
weight. Such a cost or weight function depends on t, K^, At and number 
of missions N . It is possible to select these design variables by 
minimizing the cost or weight function subject to appropriate reliability 
bounds. The effect of NDI is indirectly related to initial flaw distri- 
bution. Additional NDI effects such as the rejection of structures are 

* 

not considered in the analysis. However, they can be included as cost 
units related to the probability of failure. A numerical example is 
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illustrated in the next section to illustrate the developments of the 
report. 


Numerical Example and the Computer Program 

For the numerical example, it is assumed that the ’Johnson S. ’ 

b 

distribution for the initial crack depth is such that the minimum crack 
depth is zero and the maximum crack depth is 0.1 inch. Different possible 

ratios (a/c) are considered. Paris’ equation for crack growth is assumed 

-18 2 
with c = 0.847 x 10 and exponent equal to 3.0. The variation of v 

with (a/c) is approximated by a quadratic relation. 

The primary objective of reusing the SRM case is to reduce the 

cost of operation of the shuttle. However, as the number of uses (or 

cycles) is increased probability of failure increases because of large 

crack depths associated with more use. The probability of failure also 

increases with higher proof factors because of higher stresses. Thus, 

smaller number of cycles and small proof factors, result in higher 

reliability. However, small number of cycles increase the cost of the SRM 

case because it has to be replaced after relatively smaller number of uses. 

Then the total cost function consists of (a) the cost due to number of 

uses and proof factor and (b) the expected cost of failure, i.e., 

C total * C < S - KP > + Vf 

In the equationC(N,K^) is the cost due to number of uses N, and 
proof factor K^. The cost of failure of SRM case is- denoced by C ^ and 
the probability of failure by P^. The cost C(N,K^) can be expressed as 

1 


10 



C(N,K P ) = CjN a + c 2 (K P ) b 


It is to be noted that the expression is only for the purpose of 
illustration in this report and can be changed to reflect the figures 
more accurately. 

The power 'a 1 is negative to reflect the fact that the effective 
investment cost is lower if more number of uses can be obtained from the 
same vehicle. Similarly the power , b l is also negative'. This is to 
reflect the fact that the capability of vehicle to withstand higher proof 
load usually indicates larger available margin of safety and increased 
confidence in the success of the next mission. This also includes intan- 
gible cost due to confidence. It is to be noted that t and'^t 1 are not 
varried in the numerical example. Therefore* there is no cost associated 
directly with t or At. 

Initial thickness of the case is assumed to be 0.686 inch and it 
is assumed that 1% of the thickness is reduced after each use. The flight 


loading is assumed to be 936 psi. For the purposes of the illustrative 
example, the problem posed is to select the number of use cycles and 
proof factor for minimum expected cost. A reliability restraint can be 
imposed. However, the numerical example has not been considered such a 
restraint. Arbitrarily, the following values have been used for C^, 
and C^, c^ = c^ = 1000.00 units, c^ = 180 units, a = -0.3 and b = 4.0 
have been used. 

The general procedure can be summarized in the following steps. 

A computer program has been written to carry out the needed computations, 
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1. Obtain the parameters of the Johnson S, distribution for the 

b 

initial flaw size. 

2. Obtain the stress in the membrane from the known geometry of 
the case and wall thickness 


In the equation is the proof stress factor. During flight 
is equal to'one. Pressure P is the MEOP pressure and R 
is the radius of the SRM case equal to 72.5 inches. 

3. Obtain the new CDF and density function for the ctack depth 

p 

after the proof test. A value of K close 1.0 is assumed to 
start the calculation. 

4. Obtain the new CDF for the crack depth during the flight 
following the proof test. 

5. Estimate the probability of failure. 

6. Compute the cost function parameters. 

7. Obtain the new CDF after the material removal. 

8. Repeat steps 2 to 7 for the new thickness and the next mission 
until the total number of missions are complete. 

9. Change 1 , t y , N and repeat the calculations as necessary. 

10. Select the design variables for the minimum value of the ob- 
jective function subject to reliability constraints. 

A computer program has been written to carry out these steps. Only 
N is varied in step number 9. The program is listed in Appendix II. 
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Figure 1 illustrates the variation of cost with number of cycles 
and proof factor in the range 1.02 to 1.20. From the assumed arbitrary 
cost figures minimum expected cost occurs for 16 cycles and proof factor 
of 1.12. The corresponding reliability is only 0.9. Lower proof factor 
need to be used for higher reliability. In the numerical example pre- 
sented in this report, t and At have not been varied. 


13 
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Conclusions and Recommendations 

This report has demonstrated that the reliability analysis based on 
probabilistic fracture mechanics can be used to optimize the selection of 
the design variables of the SRM case. In particular, basic design 
variables such as the thickness and projected design life as well as the 
fracture control variables such as the proof factor and material erosion 
can be included. Accuracy in estimation of the initial flaw size distri- 
bution is reflected in the assessment of the risks involved in the design. 
By knowing the risks involved in the design, weight and cost can be 
reduced from those obtained by deterministic analysis and use of arbitrary 
safety margins. 

This report is only a first step in the development of procedures 
based probabilistic fracture mechanics. Additional work that is necessary 
can be listed as follows: 

1. A more accurate analysis can be obtained by considering the' 
joint distribution for the crack depth and crack length along 
the surface. 

2. Accurate methods of estimating the probability distribu- 
tion for the initial flaw size distribution should be developed. 

3. In particular, effects of slap down and time dependent crack 
growth including stress corrosion should be considered in the 
SRM analysis. 

4. Uncertainties in external loads and material properties should 
be considered. 
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5. Accuracy of the different models for crack growth (in the point 
of view of probabilistic fracture mechanics) should be evaluated. 

6. Alternate fracture control plans and more accurate stress inten- 
sity measures based on cylindrical geometry can be considered. 

7. Cost of NDI efforts in relation to the cost that will be 
incurred by additional safety factor should be evaluated in 
the point of view of improved reliability. 
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APPENDIX I 


Estimation of the new CDF of crack depth after use from 
knowledge of the old CDF and probability density before 



Crack Growth Rate 


The rate at which the crack depth increases is given by Paris 
equation as follows. 

= C(AK) n =_ 0.847 (dK) n x 10" 16 

For subsequent convenience in algebra, the value of f n l is taken to be 

3.0. The suggested value from current state of art is 2.48 (C = 0.847 
- 18 

x 10 ). Now substituting for AK 


da 

dN 


^ = 0.847 


C 5 + C 2 ( t } + C 3 ( t )2 


% 


Simplifying this further. 


-- = c 

dN 6 


C 5 + C 2 ( t ) + C 3 ( t )2 


where 


C 6 = 0.847 X (C 4 ) 3 X lo' 18 


da 

Integration of (— ) 
dN 


da 

Separating the variables a and N in it follows that 


dN = -zr- 

C 6 


c 5 + 0 2 (|) ■+ C 3 (|) 2 


1.5 

da 


17 
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Integrating both sides between state (1) and state (2) 


N, 


W - c7 / 

N, 6 i. 


S + Vt> + C 3 ( ! )2 


(a) 


1.5 


1.5 


da 


In order to evaluate the integral on the right hand side, it is found 

necessary to expand the numerator of the integrant binomially * 

Now consider the numerator of the integrand with C<- - 1. Neglecting 

a 3 

terms of higher order than (— ) , it follows that 


1 + vt> + vt> 2 


1.5 


= 1.0 + 1.5 C 2 (|) + I 1.5 C 3 + 1.5(0.25) | (|) 2 
+ I 1.5(0. 5) C 2 C 3 - 0.25(0.5) 2 C 2 3 J (|) 3 


Letting 


P 1 = c 1 ' 5C 2 


P = -4 1.5 C„ + 1.5(0.25) C_ 2 j 

2 c Z ( 3 Z \ 


and 


P„ = 


1.5(0. 5) 


C 2 C 3 


- (0.25) C, 


it follows that 
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1 + C 2 (f) + C 3 (f) 2 


2 3 

= 1.0 + P 1 a + P 2 a + P 3 a . 


Substituting in the integral 


[H] - 5 s - 

N l . O 


1 r \-°* 5 » _Zi / x 0 ' 5 . P 2 , . 
0.5 ^ + 0.5 ^ + 1.5 ^ 


~l 


P~ 2.5 

23 Ca) 




Solution of as a function of a^ 


Substituting the limits 

C 6 C-, - N X ) - - 2 (a 2 ,- 0 - 5 ♦ 2 P l( a 2 ) 0 ' 5 + § V/'* + f W’ 


+ 2 C a 1 )-°- 5 - 2 PjCap 0 - 5 ' - § P^) 1 ’ 5 - § P^) 2 ' 5 


Rearranging and neglecting terms of order higher than three, it reduces 
to the following 


(a^ 3 + P (a^)^ + q a L + r = 0 


where 


P = 


1.0 


r» PF i 

^ 3 12 5 r 3 J 


• [4 p --pi 

L 1 3 2 J 
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- 1.0 


q = 


fipp - — pi 
3 12 5 r 3 


[8 P L + (cp 2 ] 


and 

r = - — . 

r® p p - » p i 

.3 12 5 3 

Now, the three roots of this cubic equation, (a^) 1 are given by the 
following [CRC tables 17th edition, P. 105] 


where 



P 

3 


P 

3 


— 1 2 
a = i (3q - P Z ) 


b = 27(2P 3 - 9Pq + 27r) 
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Trans formation 


Probability density of a^ is given by 

y* 2 > 


CDF of a.^ *- s then 


W 


I £ a < a 2> ^2 ' I f a, < a l> da l 


r i a i^ a 2^ 

-[ F aM 
1 0 


where ^ a ( a \ is th e CDF D f Johnson S distribution. 

11 £> 

Now, it is needed to obtain as a function of a^, No. of 
cycles etc. This can be done by solving the polynomial equation 
obtained previously in terms of and treating a^, and N^as 
constants. The infinite degree polynomial equation is truncated at 
the 3rd degree for convenience. 

Of the three roots only one will be the real root because of 
the physical nature of the problem, say ^( a 2 ^ 

Then, substituting in the expression for the CDF of a^ 

w 

yv - J yv d <v 

z o 1 

or if the CDF of a is known, 

1 A , s 

a l (a 2 > 

F a‘ a 2> - [ F a< a l>] 

0 
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Thus F , (a 0 ) is a function of the parameters of the initial flaw 
a 2 2 

distribution i.e. X , e , T and , the proof test factor K p and 
the number of uses (N^ - N^) . 

The effect of each of these parameters can be studied by calcu- 
lating F (a ) for various cases, by means of a computer. 
a 2 1 

2 a 

Parabolic Fit to 0 (— ) 

Consider the range /<(~)<1.0. In this range it is attempted 

2 a 

to fit a parabolic curve for 0 (~) . such as follows. 


0 2 <f>-Ci +V <f) + C 3 (f)' 


In order to determine the three constants C^ 5 and three points 
are considered on the given curve. 


t - 0 


(ii) (J) = 0.5 


(if i) (“) = 1.0 


0 2 (J) = 1.0 

, 0 2 (|) = 1.5 
0 2 (|) =2.5 


Substituting the values for point (i), 


c L = 1.0 


Substituting the values for point (ii) 


1.0 + C 2 (0.5) + C 3 (0.25) = 1.5 


or 


2C 2 + c 3 = 2 -° 


22 


Substituting the values for point (iii) 

1.0 + C 2 + C 3 = 2 .5 
or 

C 2 + C 3 - 1.5 

Solving equations (2) and (3) simultaneously 

C 2 = 0.5 

and = 1.0 

Thus the chosen parabolic fit is as follows 

o 2 

0 (-) = 1.0 + 0 . 5 (— ) + (-) 

c c c 

Limits of Integration for the CDF of ’a^ 11 

By hypothesis, the initial flas *a^ * has a Johnson - distri- 
bution* Also, there is a functional relationship between the initial 
flaw size *a^ 1 and the subsequent flaw size 'a^ 1 after N cycles. This 
relationship renders 'a^ 1 a random variable because r a^ 1 is a random 
variable by hypothesis. Having known the range space of 1 the range 
space of 'a^ 1 can be derived from the functional relationship between 
'a^ 1 and 'a^ ' * Thus, if the lower limit of *a^ 1 is zero, it follows 
from the functional relationship between f a^ and f a^ * that the lower 
limit of a^ is also zero. Next, if the upper limit of *a^ * is a^ the 
corresponding upper limit for ’a^ 1 can be obtained by solving the cubic 
relation between a^ and a^, as a function of the number of cycles 
N = N 2 - Hj. 
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IAM MAIN 


74/74 0 P T=‘l 


FTN 4.6+420 


76/08/17, 




C 

c 

0 

C 


103 

C; 

102 

C 


PROGRAM MAIN ( INPUT* OUTPUT * T AP 66- INPUT i TAP E6=0UTPUT) 

8iBKIMS!illW« M ■ . . - ,* -:«¥'■' ' •. %-m . .?. 

READ (5*+ ) ALAMO A » EATA .GAMA- 
E P S - 0 • 0 

WRITE (6*100). A L AMOA , E A TA * G AM A , E PS 
ENDFILE6 " 



ltLAMBDA 
1F5.2) 

DO 50o0 1=1.10 

— AKP = 1.C+(0.02*I > 

f ?£ 1 6 J f e I ? I k ? R 00 ^' TE S r J AC T ° R M. 

FORMAT (///✓/ *5X," PROOF - TE ST 

— • DO— 40 00 — N=i * 25 — — — — — -■■ ■ ^va- . ...- r — ■ , , ,- 

Do’ , 3310 T L=1 , *2 M8ER ° F THE CYCLE CONSISTING ■■OF: ONE PROOF-TEST ANO ONE USE, 
IF(L.EQ.l) WRITE<6*104) N ■ V.vL- ...‘ •‘ :: 2 

- IF (L.EQ.2) WRIT E’(6* 105 


‘ 9 ?S SW'SSSfc i 1 ; Bi 



t g ?7 < ! J £ S g El g ^1 ^ ^ 7 


108 

C 

e - ■ 


I F ( t V E Q .-2 ) -s IG M P =S I G M P/ A KP^~ 
WRITE(6»108)SIGMP 
FORMAT (5 X* “PROOF-STRESS =“,F26 .1) 



C4= (SORT (1.2+3.147) )+SIGMP,. ; 
C6=0 . 847+(C4++3 . 0) ' 

C6*l. 0E-18+C6* 

• Pl=l. 5 + C2/C j 

P2 = 

P3 
A2 

zi 

Z2 = 




oo 


74/7/+ 0PT=1 


FTN 4.6+420 


7 6/08/17. 


Z3 = 2 . G*Pi* (SQRT ( A2 ) ) . ;; u V- k\l 

Z4 = -(2.0/SQRT(A2))**(C6*N) M 

C7=Z1+Z2+Z3+Z4 a ■ •■ "-r. % 

Z2=(8.l3*Pl*P2/3.0)-(8#QfP3/5.0) .-.v' •*.+;'• •• * '. ' • : • 

P=Z1/Z2 ■'■ ■■■.•■-':•.' ."&#. '••••''.•> -.. . i 

Zl=-< <8.0*P1)+(C7+C7) ) • • ■■ • / •;■• ■ •• ■ • y- •■ ..,.. 

R=4 0/Z2 • •. : ' +?#:■■ -V ••" y >' •• : ‘ ' ’’ V ! ^ 

ABAR=((3.0*Q)-( P+P ) ) /3 » 0 ?®li 

-z 2 = - 9 ♦ p ;q--..--. *- V .--. . . 

Zi = (B + B/4. 0) + ( ( A BAR/3. 0) + * 3 . 0) ' ! ; MV-'4 : ■ V*:“ ' >• • ■'• v •. : ;. , ..'; :* 

FOR HAT ( / * M 0 ISC REM I N ENTi-SQUA RE™«'^'rEi."5"*'5 h ^ Tr 'w 1 ' — "; ' i* — 

AA=( (-B/2.0)+CSQRT (CMP LX (Z 1, Q * 0 ) ) ) 7> - V/ 

— Q2 = A I MA G ( A AT ~ ■■■ - ' 

RR=(SQRT((Ql++2.0> + (Q2*+2.0)))++U.0/3.0iy-..,a i :. ' L' •.&*'-?'• 7>-’ ■/ '*p:> .- y' 

TH£T = ATAN (Q2/Q1 ) . •• : .y..y. - .-■.•• =: " * v yK/ . V> . ' -7 V- .:■■ ’ ‘ : .:. ;• 

AR1=RR+C0S (THET / 3. 0) ' !•■'••'• •••.'•,• ' 

■All =RR*S IN (T HE TV 3 . 0)-~— . s » 


A 1 3 - R R+-S IN t ( TH &TV~(4^3~rl 4-7-^37 0 ) i -• > ••• ; , ?:. 
Q1=.REAL(BB) - :. • . • •••>* 

Q2= AI MAG ( BB) 

RR=(SQRT({Q'1**2.0) + (Q2**2.0) >)**<!. 0/3,0) 
- THET = A TAN (Q2/Q t\— 

BR1 = RR*C0S (THET/3,0) •■V7 ; :*g 

BI 1=RR* SIN ( THET / 3. 0 ) V*?> y- 
BR2-RR*C0S C. ( THET* ( 2 . +3 . i 47 ) ) /3 » 0 
B 1 2 =R R+ S I N ( ( f H E-T-+T 2 -.-+3-.-1 

BR3 = RR+COS ( (rHET+(4.*3.147) ) /3, 0) $ S’ V 
BI 3 = RR 4 S IN ( ( THE T+ ( 4. +3 . 1 47 ) ) /3 . 0 ) : 1 
AlHAT = ARl+BRl-( P/3.0) C-'-. 

A A = REAL (A 1-HAT) ~"T 


P=Z2/Z1 • , 

Q=4.C/Z1 

Al= {.-P + SQRT < CP*P)- (4.0+Q) ) ) / 2. 0 

WRITE (6 , 114) Al -•.-• 

FORMAT (5X , "UPPER" LIMIT ' OF .Al:^. 
A 1 H A 1 = A 1 ■.!■■■■■ -v ns(i ?.%•* * ; 

IF (A1HA1, GEVALAMDA) Y (Nl»l)^i. : 
IF (A 1 H A 1 , GE-. A L AHO A > — G O-T0“3 DO 5- 




, m . pi r^i # i £ 



\S A- Lang ley, 1977 


to 

<1 


^AM MAIN 

fr~ 


74/74 0PT=1 


3 000 


3005 

116 

C 

C 


118 
3010 
4 000 
5000 
9 99 


FTN 4.6 + 420;;,,'': 


76 /Q 8 / 17 , 



THIS SECTION CALCULATES > 

ni = 15 :•* 

DO 3000 K=2,N1 
Al= A1H A1 *( K-1T/F LO ATI N 1- 1 ) 
Z1=EATA*ALAMDA/SQRT (2. 0 
Z2 = 1,Q/(A1MALAMQA-At) ) 
Z4=EATA*AL0G(A1/(ALAMDA-A1) ) 









147 ) v 

■. i , 


y-f 


’ V ; . v • M -• \-;y ■ * ■ ■ 

; ■ ■. '•<; 




I LURE = ", E 1 4. 51 ' ;; v' 


+f?******+fffff*+tff+f*++**+. 


CALL- IN T GRL M ♦ A 1HA t f F ,Y ,NB 
WRITE(6»116)Y(N1*1) T 

FORMATT5X, "PROBABILITY OF NO FAILURE = 

*****♦*♦*♦**♦***♦*♦♦*♦♦ 

•THIS 
C8 = 

CIO 

cu 

ctot=co* (cio+cii+i»o*zzi) ■ ' • a 

WRITE (6,118)CTOT ^ if ' ' -r ‘v f 

FORMAT (5X, " + **♦ + ***♦♦**♦■******** TOT AL C O S T T U N C T 1 0 H , t ’N F 1 5 ♦ 5 ) : ; ;^yvf ; , ■ ' 

CONTINUE *■ — 

continue 

CONTINUE - 
STOP 

£MD ^ 







